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Well-known Nuclear Magnetic Resonance experiments show that the time evolution according to 
(truncated) dipole-dipole interactions between n spins can be inverted by simple pulse sequences. 
Independent of n, the reversed evolution is only two times slower than the original one. Here we 
consider more general spin-spin couplings with long range. We prove that some are considerably 
more complex to invert since the number of required time steps and the slow-down of the reversed 
evolutions are necessarily of the order n. Furthermore, the spins have to be addressed separately. We 
show for which values of the coupling parameters the phase transition between simple and complex 
time-reversal schemes occurs. 



Due to wide applications of Nuclear Magnetic Reso- 
nance in medicine, chemistry, and physics much effort 
has been undertaken to understand the dynamics of nu- 
clear spins in solids and liquids . Nowadays, NMR 
plays an important role in the first experimental realiza- 
tions of quantum computation Q. 

Phcnomcnologically, the dynamical evolution of a sin- 
gle spin seems to be a relatively simple relaxation and 
dephasing process of an open quantum system, formally 
described by the Bloch equations However, the 

fact that not all dephasing and relaxation processes are 
irreversible is important for practical purposes since refo- 
cusing techniques are applied with great success J|^,||,0| ■ 

The simplest version of refocusing is possible when de- 
phasing of the (uncoupled) spins is caused by spatial in- 
homogeneities of the strength of the magnetic field. In 
this case the time evolution can be reversed by sand- 
wiching the natural dynamics by 180°-rotations of all 
spins around an axis orthogonal to the magnetic field. 
More sophisticated versions of refocusing are necessary if 
dephasing and relaxation processes are caused by spin- 
spin interactions. An important example is the so-called 
dipole-dipole coupling 0. 

We describe this interaction formally as follows. Let 
(C 2 )®" be the Hilbert space of n spins and be the 
Pauli matrix a a acting on the fc-th spin for 1 < k < n and 
a — x,y, z. Then the so-called truncated dipole-dipole 
Hamiltonian H d is given by 

k<l a 

where wu is the strength of the interaction between the 
spins k and I. 

We assume that all local unitary operations on the 
spins can be produced in arbitrarily small time (on a time 
scale governed by the strength of the coupling). This is 
often called the fast control limit. In practice this is done 
by external electro-magnetic pulses. 



Consider a pulse implementing a rotation of each spin 
around the y-axis by 90°, i. e., if the spin is in z-direction 
it is in x-direction after the rotation. Formally, one ap- 
plies a unitary transformation on (C 2 )®" of the form 

Vy := Uy ® Uy <g> • • • Cg) Uy 

with Uy(j z Uy — a x . Now we focus on the effect of the fol- 
lowing scheme: (1) implement v y by the corresponding 
pulse, (2) wait the time t, and (3) implement v y , i. e., the 
inverse of v y . The net effect of the scheme is an evolution 

Vy exp(-iH d t)v y = exp(-ivlH d v y t) , 

i.e., a dynamics according to the Hamiltonian v y H d v y . 
We choose a second rotation u x around the x-axis such 
that u x a z u. x = a y and set v x := u x ® u x (£> • ■ • (£> u x . It is 
decisive to note that 

-H d = v y HdVy + v\H d v x . 

This has the interesting consequence that for small times 
e one has 

v x exp(-iH d e)v x exp(—iH d e)v y w exp(+iH d e) 

up to terms of second order of e. This first order ap- 
proximation is usually referred to as average Hamiltonian 
theory [Q. The intuitive meaning is that the time evo- 
lution according to H d can be approximative^ inverted 
by dividing the time t into small intervals of length e 
and interspersing the natural time evolution by the above 
pulses. The inverted dynamics is two times slower than 
the natural one. We express it by saying that the time 
overhead of inverting the Hamiltonian H d is 2. The task 
of time-reversal is a special case of the general problem 
of simulating Hamiltonians that has recently been con- 
sidered by several authors p|-|l9|]. 

Note that the decoupling and inversion schemes for 
the dipole-dipole Hamiltonian is simple with respect to 
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three aspects: (a) it does not use selective pulses, i.e., 
in each time step the same unitary transformation is ap- 
plied to all spins, (b) The number of steps of this inver- 
sion scheme does not increase with n since it is always 2, 
and (3) the time overhead does not increase with n since 
it is also 2. 

The purpose of this article is to show that there are 
interactions that are considerably more complex to in- 
vert with respect to all three criteria. The most general 
spin-spin interaction on n spins is given by 

k<l a, 13 

where J is chosen to be a real symmetric 3n x 3n-matrix 
with zeros for k = I. Note that the symmetry of the cou- 
pling matrix J does not imply any physical symmetry 
of the interaction. It is a consequence of our redundant 
notation that turns out to be very useful. The coupling 
matrix J consists of 3 x 3-blocks. The 3 x 3-matrix Ju 
given by the block at position (fc, I) describes the cou- 
pling between the spins k and I, We have Jik = J^, i.e. 
the matrix describing the coupling between the spins I 
and k is just the transpose of the matrix describing the 
coupling between the spins k and I. The blocks on the 
diagonal are zero matrices. 

Let C := SU(2) <g> SU(2) <g> • • • ® SU{2) be the con- 
trol group of local operations on the spins. As already 
mentioned all operations in C are assumed to be im- 
plemented arbitrarily fast. The most general refocusing 
scheme (also called time-reversal or inversion) based on 
the first order approximation above is given by rotations 
vi, v%, . . . , Vn £ C, and (relative) times ti, t%, . . . , tjv such 
that the average Hamiltonian 

//: Y^'.rlljr, 

3 

is equal to —Hj. Here N is the number of time steps and 
t := tj is the time overhead. The time overhead 

gives the slow-down of the inverted Hamiltonian. In the 
following we consider both the number of time steps and 
the time overhead as complexity measures. In contrast, 
if the sequences are chosen in such a way that the av- 
erage Hamiltonian is zero one has a decoupling scheme. 
Each inversion scheme with N steps defines a decoupling 
scheme with N + 1 steps in a straightforward way (note 
that the notion of time overhead does not make sense for 
decoupling). Conversely, each decoupling scheme can be 
converted to an inversion scheme as follows. The equa- 
tion 

JV 
3=0 

implies 

JV 
2=1 



by elementary calculation. Therefore we shall restrict our 
attention to inversion. 

The condition on inversion schemes can be rewritten 
by expressing the action of the control operations on the 
coupling matrix J as follows. Note that any unitary op- 
eration u £ SU(2) corresponds to a rotation on the 3- 
dimensional Bloch sphere via the relation 

a a 

where the vector c = (c Xl c y ,c z ) is obtained by applying 
a rotation U £ SO(3) on the vector c — {c x ,c y ,c z ). It 
is straightforward to verify that conjugation of Hj by 
v := w 1 ' ® (g) • • • (g) u' n ) corresponds to conjugation 
of J by a block diagonal matrix of the form 

n 

V := C/ (1) © U {2) © ■ • ■ © U {n) G 50(3) . 

fe=i 

The condition for time inversion is hence given by 

-J = J2 t i V i JV ?> (!) 

3 

where the orthogonal matrix Vj corresponds to the uni- 
tary transformation Vj for j — 1, . . . , N. Time reversal 
and decoupling schemes that apply for general coupling 
J have been presented in . These schemes can even 
be used when the coupling J is unknown. Here we fo- 
cus on optimality criteria for schemes referring to spe- 
cific couplings and show that there are interactions that 
cannot be inverted significantly more efficiently than an 
unknown interaction. 

We restrict our attention to interactions with an addi- 
tional symmetry, namely a Hamiltonian of the following 
form 

k<l a,0 

The matrix W := (uiki) is a real symmetric n x n-matrix 
with zeros on the diagonal. It describes the coupling 
strengths and the signs of the interactions between all 
spins. The matrix A = (A a p) is a real symmetric 3x3- 
matrix characterizing the type of the coupling. This 
means that all spins interact with each other via the same 
interaction and that only the coupling strength and the 
sign varies. It is important that in this special case the 
coupling matrix J can be expressed as a tensor product 
of W and A, i.e. J — W <g> A. 

The matrix W characterizes the coupling topology of 
the spins. In the language of graph theory the matrix 
W is the adjacency matrix of a weighted graph. There- 
fore many results on graph spectra (eigenvalues of the 
adjacency matrices) can be used to derive lower and up- 
per bounds on the number of time steps and the time 
overhead. 

To discuss the complexity aspects we distinguish be- 
tween the following three cases: 
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1. A is traceless. 

All spins can be subjected to the same transforma- 
tions in each time step, the number of time steps 
and the time overhead are at most 2. 

2. A has negative and positive eigenvalues but 
tr(A) ^ 0. 

The spins have to be addressed separately, the num- 
ber of time steps necessarily grows for increasing n. 
But the time overhead does not depend on n. It de- 
pends only on the eigenvalues of A. 

3. A is either positive or negative semidefinite, i.e., 
the non-zero eigenvalues have the same sign. 

Then the spins have to be addressed separately, the 
number of time steps is at least n — 1, and the time 
overhead is also at least n — 1 . 

To prove these statements we assume w.l.o.g. that the 
interaction between all pairs is of the form 

a x cr x ®a x + a y a y (g) a y + a z a z <g a z , 

where a x ,a y , a z are the eigenvalues of A. This can always 
be achieved by rotating the reference frame (l2) . 

Case 1. Assume A to be traceless. Let 5 be a rotation 
in the Bloch sphere that realizes the cyclic permutation 
of the axis according to x — > y — > z — > x. Then we have 

2 

^V j JV jT = 0, 

j=o 

where V := 5 © 5 © • ■ • © 5 and V j is the jth power 
of V. Here 5 is a rotation in the Bloch sphere that 
realizes the cyclic permutation of the axis according to 
x — ► y — > z — » x. The jth power of V is denoted by 
V' } . Hence the Hamiltonian is inverted by a sequence of 
length 2. Due to the equation 

-J = VJV T + v 2 jv 2T 
the time overhead of this inversion scheme is 2. 

Case 2. The fact that the spins have to be addressed 
separately has also been noted in . It can be seen as 
follows. If all spins are subjected to the same transfor- 
mation each 3 x 3-block of J is conjugated by an element 
of 50(3). Such a conjugation preserves the trace of each 
block matrix. This trace is even preserved under posi- 
tive linear combination, i.e. the trace of each resulting 
3x3 matrix has the same sign as the trace of the original 
matrix. Therefore a sequence applying the same trans- 
formation to all spins can never change the sign of the 
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trace in any block. But this is required for obtaining the 
inverse entries. 

We show that the number of time steps cannot be 
constant, i. e. independent of n. Assume w.l.o.g. that 
tr(A) > 0. Choose a partition of 50(3) = 5iU5 2 U. . .U5 P 
into equivalence classes such that for any pair 0, G 5j 
(i = 1, . . . ,p) implies tr(OAO T ) > (p depends on A). 
This is always possible since the partition can be chosen 
in such a way that all orthogonal transformations in the 
same class are sufficiently close together. Then OAO T 
almost corresponds to a conjugation of A by an orthog- 
onal matrix and its trace cannot deviate too much from 
tr(A) due to the continuity of the trace. 

Let Oj k ^ £ 50(3) be the operation performed on the 
kth spin in time step j (described as action on the cou- 
pling matrix). In each time step at least \n/p] spins are 
transformed by operations of the same class 1 . After N 
time steps at least \n/p ] spins have been transformed 
by operations of the same class in each step. Let k and I 
be such a pair that has been transformed by operations 
of the same equivalence class in each time step. Then the 
trace of the average coupling between the spins k and I 
after the sequence is given by 

*k5> ?° ao ? t )>®- ( 2 ) 

3 

If the sequence is an inversion scheme then we must have 
^ tj 0fAof = -A. 

3 

But this condition contradicts the inequality (|^). There- 
fore such a pair k and I must not exist. This gives the 
lower bound N > log nj logp on the number of time steps 
N. 

We show that the time overhead does not depend on 
n. We assume w.l.o.g. a x > > a z for the eigenval- 
ues of A. First we describe a partial decoupling scheme 
selecting for instance the a z © a z coupling terms while 
switching off the xx and yy terms. Following jl4| such a 
decoupling can be achieved by certain sequences of local 
conjugations by the unitary ia z . Each time step of this 
scheme is described by a column of a Hadamard matrix. 
The entries determine which spins are conjugated by ia z 
transformations. The idea is that the xx and yy terms 
acquire in exactly half of the time steps a minus sign. 
Note that this scheme does not weaken the zz coupling. 
In other words, although this scheme requires a number 
of time steps of the order n the time overhead for the sim- 
ulation of the interaction ^ Wkia z a z k ^ 'ci'' by the original 
one is 1. 
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There are two cases depending on whether a y > or 
a y < 0. Consider the case a y > w.l.o.g. First switch 
off the xx and yy interactions. Apply on each spin a 
conjugation by u y (defined as in the beginning). This 
simulates the interaction 

k<l 

By applying this interaction for the (relative) time 
—dx/iz we have reversed the ccir-components in the orig- 
inal Hamiltonian. Since a y is not negative we can reverse 
similarly the yy-components with time overhead —a y /a z . 
To invert the zz-components we use the ra-components. 
The time overhead is —a z /a x . In summary we see that 
the time overhead is independent of n. It only depends 
on the eigenvalues of A. 

Case 3. Let t 1; f 2 , ■ ■ ■ , i/v be the (relative) times and 
V\, V%, . . . , Vn with Vj 6 ©fc = i SO (3) be the operations 
of an inversion scheme, i. e., 

*j Vi {W ®A)Vf = -W ®A. (3) 

3 

First note that an inversion scheme for a given Hamil- 
tonian Hj can also be applied to a rescaled Hamiltonian 
obtained by changing the strength and the sign of an ar- 
bitrary spin pair interaction. If all the coefficients wu 
are non-vanishing for k ^ I then we can equivalently de- 
scribe a time inversions scheme for the coupling matrix 
J := W <E> A where all non-diagonal entries of W are 1 
and the diagonal consist of zeros. Hence the fact that in- 
teractions decrease with distance of the spins is irrelevant 
for considerations of the complexity of time inversion (as 
long as the interaction cannot be neglected). Therefore 
we assume w.l.o.g. all non-diagonal entries of W to be 1. 

For each time step let Vj be the block diagonal matrix 

Vj = Uf ] © • ■ • © <7j n) . Set R := £f =1 tjVj(l © A)Vj . 
We add R to both sides of eq. @ and obtain 

N 

^2 tjVj((W + 1) © A)V J T = —W <S> A + R. (4) 

The rank of the matrix (W + 1) is 1 since all its en- 
tries are 1. Consequently, the rank of the left hand side 
of equation (Q) is at most N r, where r is the rank of A. 
Note that R is a semi-positive matrix. Therefore the rank 
of the right hand side is at least the number of positive 
eigenvalues of — W<E>A. The matrix —W has n—1 positive 
eigenvalues (by adding the identity to W one easily ver- 
ifies that W has the eigenvalues n — 1, — 1, — 1,...,— 1). 
Hence — W © A has [n — l)r positive eigenvalues. We 
conclude that the rank of the right hand side is at least 
(n — l)r. This proves the statement N > n — 1. 

We show now that the time overhead also grows with 
n. Let A m in and A max be the minimal and maximal eigen- 
values of J, respectively. Note that A m ; n is negative since 



J is traceless. Let t\,t^ . be the (relative) lengths 
of the time steps. Then we have 

^ ^ tj Amax 111 ^min ; 
3 

since the maximal eigenvalue of the sum of matrices in 
eq. (|l|) is at most the sum of their maximal eigenvalues 
p0[ and — A m i n is the maximal eigenvalue of — J. We con- 
clude that the time overhead r := J^j tj f° r the inverse 
evolution is at least 

T ^ A max/ Amin ■ 

Since the maximal and minimal eigenvalues of W are 
(n — 1) and —1, respectively, the maximal and minimal 
eigenvalues of J are (n — l)a max and — a max where a max 
is the maximal eigenvalue of A. We conclude that the 
minimal time overhead is at least n — l.D 

The truncated dipole-dipole coupling is an example 
for case 1. The corresponding matrix A is given by 
A := diag (1, 1, —2). An important example for case 3 is 
the strong scalar coupling where A is the identity matrix. 
By combining both types one can easily find examples for 
case 2 

Besides the practical applications of refocusing the 
lower bounds on the inversion complexity may be seen 
as a small contribution to the old debate on the origin 
of phenomcnological irreversibility: Whereas the consid- 
ered dynamical evolutions are strictly reversible the com- 
plexity of a real implementation of the time inversion is 
growing with the number of interacting particles. Note 
that the relevance of complexity theory for the defini- 
tion of physical entropy has been advocated for by sev- 
eral authors |p2]-p4]|. Of course our complexity bounds 
are restricted to the assumptions of average Hamiltonian 
theory. Inversion schemes that are not based on this ap- 
proach and include higher order terms of the Magnus 
expansion may be less complex. However, our results 
present an example how to develop a complexity theory 
of time reversal in physics. 
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